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ABSTRACT 

We report new precision measurements of the properties of our Galaxy's supermassive black hole. 
Based on astrometric (1995-2007) and radial velocity (2000-2007) measurements from the W. M. Keck 
10-meter telescopes, a fully unconstrained Keplerian orbit for the short period star SO-2 provides values 
for the distance (Rq) of 8.0 ± 0.6 kpc, the enclosed mass (Mb,i) of 4.1 ± 0.6 X IO^Mq, and the black 
hole's radial velocity, which is consistent with zero with 30 km/s uncertainty. If the black hole is 
assumed to be at rest with respect to the Galaxy (e.g., has no massive companion to induce motion), 
we can further constrain the fit and obtain Rq — 8.4 ± 0.4 kpc and Mb/j — 4.5 ± 0.4 x lO^M©. More 
complex models constrain the extended dark mass distribution to be less than 3-4 x IQ^Mq within 0.01 
pc, ~100x higher than predictions from stellar and stellar remnant models. For all models, we identify 
transient astrometric shifts from source confusion (up to 5x the astrometric error) and the assumptions 
regarding the black hole's radial motion as previously unrecognized limitations on orbital accuracy 
and the usefulness of fainter stars. Future astrometric and RV observations will remedy these effects. 
Our estimates of Rq and the Galaxy's local rotation speed, which it is derived from combining Rq with 
the apparent proper motion of Sgr A*, (6*0 = 229 ± 18 km s~^), are compatible with measurements 
made using other methods. The increased black hole mass found in this study, compared to that 
determined using projected mass estimators, implies a longer period for the innermost stable orbit, 
longer resonant relaxation timescales for stars in the vicinity of the black hole and a better agreement 
with the Mbft,-(T relation. 

Subject headings: black hole physics - Galaxyicenter — Galaxyikinematics and dynamics — in- 
frared:stars - techniques: high angular resolution 



1. INTRODUCTION 

Ever since the discovery of fast moving (v > 1000 km 
s~^) stars within 0."3 (0.01 pc) of our Galaxy's central 
supermassive black hole (Eckart & Genzel 1997; Chez et 
al. 1998), the prospect of using stellar orbits to make pre- 
cision measurements of the black hole's mass (M^/j) and 
kinematics, the distance to the Galactic center (Rq) and, 
more ambitiously, to measure post-Newtonian effects has 
been anticipated (Jaroszynski 1998, 1999; Salim & Gould 
; 1999; Fragile & Mathews 2000; Rubilar & Eckart 2001; 

Weinberg, Milosavlejic & Chez 2005; Zucker & Alexan- 
' der 2007; Kraniotis 2007; Will 2008). An accurate mea- 
' surement of the Galaxy's central black hole mass is useful 
for putting the Milky Way in context with other galax- 
ies through the apparent relationship between the mass 
of the central black hole and the velocity dispersion, a, 
of the host galaxy (e.g., Ferrarese & Merrit 2000; Geb- 
hardt et al. 2000; Tremaine et al. 2002). It can also 
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be used as a test of this scaling, as the Milky Way has 
the most convincing case for a supermassive black hole 
of any galaxy used to define this relationship. Accurate 
estimates of Rq impact a wide range of issues associated 
with the mass and structure of the Milky Way, including 
possible constraints on the shape of the dark matter halo 
and the possibility that the Milky Way is a lopsided spi- 
ral (e.g., Reid 1993; OUing & Merrifield 2000; Majewski 
et al. 2006). Furthermore, if measured with sufficient ac- 
curacy (~1%), the distance to the Galactic center could 
influence the calibration of standard candles, such as RR 
Lyrae stars, Cepheid variables and giants, used in es- 
tablishing the extragalactic distance scale. In addition 
to estimates of Mi^h and Rq, precision measurements of 
stellar kinematics offer the exciting possibility of detect- 
ing deviations from a Keplerian orbit. This would allow 
an exploration of a possible cluster of stellar remnants 
surrounding the central black hole, suggested by Morris 
(1993), Miralda-Escude & Gould(2000), and Freitag et 
al. (2006). Estimates for the mass of the remnant cluster 
range from 10"' — lO^A/© within a few tenths of a parsec of 
the central black hole. Absence of such a remnant cluster 
would be interesting in view of the hypothesis that the 
inspiral of intermediate-mass black holes by dynamical 
friction could deplete any centrally concentrated cluster 
of remnants. Likewise, measurements of post-newtonian 
effects would provide a test of general relativity, and, ul- 
timately, could probe the spin of the central black hole. 

Tremendous observational progress has been made over 
the last decade towards obtaining accurate estimates of 
the orbital parameters of the fast moving stars at the 
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Galactic center. Patience alone permitted new astro- 
metric measurements that yielded the first accelerations 
(Ghez et al. 2000; Eckart et al. 2002), which suggested 
that the orbital period of the best characterized star, SO- 
2, could be as short as 15 years. The passage of more time 
then led to full astrometric orbital solutions (Schodel et 
al. 2002, 2003; Ghez et al. 2003, 2005a), which increased 
the implied dark mass densities by a factor of 10* com- 
pared to earlier velocity dispersion work and thereby so- 
lidified the case for a supermassive black hole. The ad- 
vent of adaptive optics enabled radial velocity measure- 
ments of these stars (Ghez et al. 2003), which permitted 
the first estimates of the distance to the Galactic center 
from stellar orbits (Eiscnhauer et al. 2003, 2005). 

In this paper, we present new orbital models for SO- 
2. These provide the first estimates of the distance to 
the Galactic center and limits on the extended mass dis- 
tribution based on data collected with the W. M. Keck 
telescopes. The ability to probe the properties of the 
Galaxy's central supermassive black hole has benefit- 
ted from several advancments since our previous report 
(Ghez et al. 2005). First, new astrometric and radial 
velocity measurements have been collected between 2004 
and 2007, increasing the quantity of kinematic data avail- 
able. Second, the majority of the new data was obtained 
with the laser guide star adaptive optics system at Keck, 
improving the quality of the measurements (Ghez et al. 
2005b; Hornstein et al. 2007). These new data sets are 
presented in ij2l Lastly, new data analysis has improved 
our ability to extract radial velocity estimates from past 
spectroscopic measurements, allowing us to extend the 
radial velocity curve back in time by two years, as de- 
scribed in ^ The orbital analysis, described in iJH iden- 
tifies several sources of previously unrecognized biases 
and the implications of our results are discussed in 

2. OBSERVATIONS & DATA SETS 

2.1. High Angular Resolution Imaging: Speckle and 
Adaptive Optics 

For the first eleven years of this experiment (1995- 
2005), the proper motions of stars orbiting the center of 
our Galaxy were obtained from K[2.2 /im]-band speckle 
observations of the central stellar cluster with the W. M. 
Keck I 10-meter telescope and its facility near-infrared 
camera, NIRC (Matthews & Soifer 1994; Matthews et 
al. 1996). A total of 27 epochs of speckle observations 
are included in the analysis conducted in this paper, of 
which 22 have been reported in earlier papers by our 
group (Ghez et al. 1998, 2000, 2005a). Five new speckle 
observations, between 2004 April and 2005 June, were 
conducted in a similar manner. In summary, during each 
observing run, ^10,000 short {t^xp = 0.1 sec) exposure 
frames were obtained with NIRC in its fine plate scale 
mode, which has a scale of 20.46 ± 0.01 mas pixel"^ (see 
Appendix B) and a corresponding field of view of 5'.'2 
X 5'.'2. Interleaved with these observations were similar 
sequences on a dark patch of sky. From these data, we 
produce images that are diffraction-limited {6 = 0'.'05) 
and have Strehl ratios of ^0.05. 

With the advent of laser guide star adaptive optics 
(LGSAO) in 2004 on the 10 m W. M. Keck II telescope 
(Wizinowich et al. 2006; van Dam et al. 2006), we 
have made measurements of the Galaxy's central stel- 
lar cluster with much higher Strehl ratios (Ghez et al. 



2005b). Between 2004 and 2007, nine LGSAO data sets 
were taken using the W. M. Keck II facility near-infrared 
camera, NIRC2 (P.I. K. Matthews), which has an average 
plate scale of 9.963 ± 0.006 mas pixel"^ (see Appendix 
C) and a field of view of 10'.'2 x 10'.'2. AU but one of the 
observations were obtained through a K' (Ao=2.12 ^m, 
AA=0.35 /zm) band-pass filter, with the remaining one 
obtained through narrow band filters (CO: Ao — 2.278 
/xm, AA = 0.048 and Kcont: Ao = 2.27 /im, AA = 
0.030 /im). During these observations, the laser guide 
star's position was fixed to the center of the camera's 
field of view and therefore moved when the telescope was 
dithered. While the laser guide star is used to correct 
most of the important atmospheric aberrations, it does 
not provide information on the tip-tilt term, which, for 
all our LGSAO observations (imaging and spectroscopy), 
was obtained from visible observations of USNO 0600- 
28577051 (R = 13.7 mag and ArsgrA* = 19")- Details 
of the observing setup for 2004 July 26, 2005 June 30, 
and 2005 July 31 are described in detail in Ghez et al. 
(2005b), Lu et al. (2008), and Hornstein et al. (2007), 
respectively. While each of these early LGSAO observa- 
tions had a slightly different setup and dither pattern, the 
more recent, deeper, LGSAO measurements (2006-2007) 
were obtained with nearly identical setups. Specifically, 
we used a 20 position dither pattern with randomly dis- 
tributed (but repeatable) positions in a 0'.'7 x 0'.'7 box 
and an initial position that placed IRS 16NE on pixel 
(229, 720) at a sky PA set to 0.0. This setup keeps 
the brightest star in the region, IRS 7 (K=6.4), off the 
field of view at all times. At each position, three expo- 
sures, each composed of 10 coadded 2.8 sec integrations, 
were obtained; the integration time was set with the aim 
of keeping the detector's response linear beyond the full 
width at half maximum (FWHM) point for the brightest 
(K=9.0) star in the field of view; the number of images 
per position was chosen to provide the minimum elapsed 
time needed to allow the LGSAO system's optimization 
algorithm to converge (~3 min.) before dithering. Table 
[T] summarizes all the new imaging data sets. 

2.2. Adaptive Optics Spectroscopy 

To monitor the line-of-sight motions of stars orbiting 
the center of our Galaxy between the years 2000 and 
2007, high angular resolution spectroscopic observations 
of stars in the Sgr A* stellar cluster were taken with both 
the natural guide star adaptive optics (NGSAO; Wiz- 
inowich et al. 2000) system (2000-2004) and the LGSAO 
system (2005-2007) on the W. M. Keck II 10 m telescope. 
The NGSAO atmospheric corrections and the LGSAO 
tip-tilt corrections were made on the basis of visible ob- 
servations of USNO 0600-28579500 (R = 13.2 mag and 
Ar - 30") and USNO 0600-28577051 (R = 13.7 mag 
and Ar ^ 19"), respectively. While the angular resolu- 
tion of the NGSAO spectra was typically 2-3 times the 
diffraction limit {Odiff = 54 mas), a point spread func- 
tion (PSF) FWHM of 70 mas at 2 ^m was achieved 
for the LGSAO long exposure spectra. 

Three different spectrometers have been used over the 
course of this study. Our earliest measurements were 
obtained in 2000 June with NIRSPEC (McLean et al. 
1998, 2000) in its low resolution slit spectrometer mode 
(R ~ 2600). It was not originally designed to go behind 
the adaptive optics system and therefore had inefficient 
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TABLE 1 

Summary of New Keck Imaging Observations 
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FWHM 


Strelil 

O LI Ulii 


l\Tl 1 VYI I'^OT' 










Obtained 


Used 


(sec) 


(mas) 




or otars 


mag 


mas 


2004 April 29-30 


Speckle 


20140 


1444 


1 X 0.137 


63 


0.09 


163 


15.9 


0.9 


2004 July 25-26 


Speckle 


14440 


2156 


1 X 0.137 


61 


0.07 


165 


15.9 


0.9 


2004 August 29 


Speckle 


3040 


1300 


1 X 0.137 


60 


0.08 


138 


15.7 


1.0 


2005 April 24-25 


Speckle 


15770 


1677 


1 X 0.137 


60 


0.07 


143 


15.6 


0.9 


2005 July 26-27 


Speckle 


14820 


1825 


1 X 0.137 


62 


0.05 


116 


15.5 


1.2 


2004 July 26 


LGSAO{l) 


12 


12 


50 X 0.181 


60 


0.31 


233 


16.0 


0.3 


2005 June 30 


LGSAO{2) 


10 


10 


5 X 7.2/11.9*' 


61 


0.32 


269 


16.4 


1.3 


2005 July 30-31 


LGSAO(3) 


66 


32 


10 X 2.8 


62 


0.34 


565 


19.0 


0.19 


2006 May 3 


LGSAO{4) 


153 


107 


10 X 2.8 


58 


0.30 


562 


19.2 


0.16 


2006 June 20-21 


LGSAO(4) 


295 


152 


10 X 2.8 


57 


0.33 


580 


19.1 


0.10 


2006 July 17 


LGSAO(4) 


70 


64 


10 X 2.8 


59 


0.31 


574 


19.2 


0.19 


2007 May 17 


LGSAO(4) 


103 


77 


10 X 2.8 


58 


0.35 


566 


19.1 


0.21 


2007 May 20 


LGSAO(l) 


20 


12 


10 X 2.8 


77 


0.20 


394 


17.8 


0.28 


2007 Aug 10, 12 


LGSAO{4) 


142 


79 


10 X 2.8 


57 


0.32 


553 


19.1 


0.20 



^ For the LGSAO data sets, the number in parentheses denotes the observational setup used (e.g., dither pattern and camera orientation; 
see §2.1 for details).^ The number of stars detected within 3 arcsec of SgrA*.*^ Kn?n is the magnitude at which the cummulative 
distribtuion function of the observed K magnitudes reaches 90% of the total sample size.^ For this analysis only stars in 4 or more 
epochs arc considered to eliminate any spurious source detections.^ The average positional uncertainty due to centroiding in each epoch 
is estimated from a set of 25 stars detected in all epcohs and brighter than K~13 mag.* Half of the images were taken using a narrow 
band CO filter, with the shorter exposure time, and the other half using a narrow band Kcont filter, with the longer exposure time. 



throughput in its AO mode; it was, however, the only 
spectrometer available behind the AO system in 2000. 
While the resulting low signal to noise data set yielded 
no line detections in the initial analysis of SO-2 (Gezari 
et al. 2002) , we now have the advantage of knowing what 
type of lines are present in the spectra and have there- 
fore included this data set in our analysis by retroactively 
identifying the Br7 line, which is used to measure radial 
velocities (see §3.2) Between 2002 and 2005, NIRC2 (P.I. 
K. Matthews) was used in its spectroscopic R ~ 4000 
mode, which is generated with a 20 mas pixel scale, a 
medium-resolution grism and a 2 pixel slit. In 2002, 
this produced the first line detection in SO-2 (Ghez et 
al. 2003) and, since then, three new NIRC2 measure- 
ments (2 with NGSAO and 1 with LGSAO) have been 
obtained. Since 2005, OSIRIS, which is an integral field 
spectrograph with a 2 fim spectral resolution of ~ 3600 
(Larkin et al. 2006), has been used. The field of view 
of this spectrograph depends on the pixel scale and fil- 
ter. Most of the OSIRIS observations were taken using 
the 35 mas pixel scale and the narrow band filter Kn3 
(2.121 to 2.229 /im; includes Br7), which results in a 
field of view of l."12 x 2. "24, and were centered on SO-2. 
All of the OSIRIS observations were obtained with the 
LGSAO system. Table [2] summarizes the details of the 
10 new spectroscopic measurements of SO-2 that were 
made between the years 2003 and 2007 (see Gezari et 
al 2002 & Ghez et al. 2003 for details of the 2000-2002 
measurements) . 

3. DATA EXTRACTION 

3.1. Image Analysis & Astrometry 

The individual speckle and adaptive optics data frames 
are processed in two steps to create a final average image 
for each of the 34 imaging observing runs. First, each 
frame is sky-subtracted, fiat fielded, bad-pixel-corrected, 
corrected for distortion effects and, in the case of 
the speckle data, resampled by a factor of two; the 
distortion correction applied to the NIRC2/LGSA0 
data is from the NIRC2 pre-ship review results 



(http : //www2 . keck . hawaii . ed.u/inst/nirc2/preship_testing . p 

and those applied to the speckle data sets are the com- 
bined transformations given in Ghez et al. (1998) and 
Lu et al. (2008). The frames are then registered on 
the basis of the position of IRS 16C, for the speckle 
images, and a crosss-correlation of the entire image, for 
the LGSAO image, and combined. For the adaptive 
optics data sets, the frames whose PSF has a FWHM 
< 1.25 X FWHMmin, where FWHMmin is the minimum 
observed FWHM for each epoch and which typically 
includes ~70% of the measured frames, are combined 
with a weighted average with weights set equal to their 
strehl ratios. To increase the signal to noise ratio of 
the 2005 June data set, the data taken through the 
two narrow-band filters are averaged together. For the 
speckle data set, only the best ^ 2,000 frames from each 
observing run are combined using a weighted "Shift- 
and-Add" technique described by Hornstein (2007). 
The selected frames from each observing run (speckle 
and LGSAO) arc also divided into three independent 
subsets from which three subset images are created in 
a similar manner to the average images; these subset 
images are used to assess photometric and astrometric 
measurement uncertainties. Figure [T] shows examples 
of the final average LGSAO and speckle images. While 
all the images sets have point spread function (PSF) 
cores that are nearly diffraction-limited {0 ~ 0'.'06 vs. 
(^diff. Urn = 0'.'05), the LGSAO images have much higher 
image quality than the speckle images, with median 
Strehl ratios of -0.3 and 0.07, for the LGSAO and 
speckle images, respectively. 

Point sources are identified and characterized in each 
of the images using the PSF fitting program StarFinder 
(Diolaiti et al. 2000) on both the average images and 
the subset images. StarFinder iteratively generates a 
PSF based on user selected point sources^ in the image 

* In this analysis, the stars that are input into the PSF construc- 
tion are IRS 16C, 16NW, and S2-17 for the speckle images and IRS 
16C, 16NW, 16NE, 16SW, 33E, 33W, 7, 29N, and GEN+2.33+4.60 
for the LGSAO images. 
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TABLE 2 

Summary of New Keck Spectroscopic Observations 



Date 


Instrument/ 


Filter: Spectral 


Fix. 


Num. Exp. 


SNR*^ 


Calibration Stars 


(UT) 


AO System 


Range (/^m) 


Scale 


X Texp (sec) 




(G2/A0) 


2003 June 08 


NIRC2'' / NGS 


K: 2.08 - 2.34 


20 


2 X 1200 


62 


HD 193193 / HD 195500 


2004 June 22 


NIRC2'' / NGS 


K': 2.00 - 2.26 


20 


16 X 1200 


23 


HD 193193 / HD 195500 


2005 May 30 


NIRC2^ / LGS 


K': 2.00 - 2.26 


20 


7 X 1200 


31 


HD 198099 / HD 195500 


2005 July 03 


OSIRIS / LGS 


Kbb: 1.97 - 2.39 


20 


7 X 900 


30 


HD 193193 / HD 195500 


2006 May 23 


OSIRIS / LGS 


Kbb: 1.97 - 2.39 


35 


4 X 900 


25 


HD 193193 / HD 195500 


2006 June 18 


OSIRIS / LGS 


Kn3: 2.121 - 2.229 


35 


10 X 900 


52 


HD 198099 / HD 195500 


2006 June 30 


OSIRIS / LGS 


Kn3: 2.121 - 2.229 


35 


9 X 900 


33 


HD 193193 / HD 195500 


2006 July 1 


OSIRIS / LGS 


Kn3: 2.121 - 2.229 


35 


9 X 900 


60 


HD 150437 / HD 155379 


2007 May 21 


OSIRIS / LGS 


Kn3: 2.121 - 2.229 


35 


3 X 900 


28 


HD 198099 / HD 195500 


2007 July 18-19 


OSIRIS/LGS 


Kn3: 2.121 - 2.229 


35 


2 X 900 


22 


HD 193193 / HD 195500 



^ For the NIRC2 data sots, the slit position angle was 259.4° (2003), 333.76° (2004), and 355.9° (2005). Tiie SNR is per 
spectral pixel and is calculated between 2.13 and 2.145 /^m. The width of a spectral pixel is roughly 2.5 and 2.53 A for 
OSIRIS and NIRC2, respectively. 
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Fig. 1. — A comparison of raw images obtained with LGSAO and speckle imaging with the Keck 10 m telescopes. The large scale image 
is an LGSAO image obtained in 2005. The inset LGSAO image (top right) and speckle image (bottom right) are centered on the black 
hole, SgrA* (marked with a cross), with a field of view of 1"0 X 1"0, also obtained in 2005. The image quality, depth, and astrometric 
precision have all been greatly improved with the advent of LGSAO. 



and identifies additional sources in the image by cross- 
correlating the resulting PSF with the image. The initial 
source list for each image is composed only of sources 
detected in the average images with correlation values 
above 0.8 and in all three subset images with correla- 
tion values above 0.6. Eleven bright (K<14 mag), non- 
variable sources establish the photometric zero points for 
each list based on measurements made by Rafelski et al. 
(2007; IRS 16C, IRS 16SW-E, S2-17, Sl-23, Sl-3, Sl-4, 
S2-22, S2-5, Sl-68, SO-13, Sl-25). As shown in Figured 
the deep LGSAO images (Kji^ ^19 mag) are three mag- 
nitudes more sensitive than the speckle images (Kji^ ~16 
mag), which results in roughly three times more sources 



being detected in the LGSAO images than the speckle 
images over a comparable region. Because of the higher 
signal to noise, as shown in Figure [31 the centroiding 
uncertainties {6X' , SY'), which are estimated from the 
RMS error of the measurements in the three subset im- 
ages, are a factor of 6 more precise for the deep LGSAO 
data sets (0.17 mas) than the speckle data sets (1.1 mas), 
for bright stars (K<13 mag); the plateau observed in the 
relative centroiding uncertainties for the brighter stars 
(K<13) in the LGSAO images is likely caused by the 
combined effects of differential tip-tilt jitter and residual 
optical distortions across the field of view. 
The sources identified each night are matched across 
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Fig. 2. — Comparison of the sensitivity of the average images from each epoch. The recent LGSAO images, with significantly longer 
on-sky integration times (ttot ~ 50 vs. 3 min) and much higher strehl ratios, are three magnitudes more sensitive than any of the speckle 
images. 



multiple epochs and their positions are transformed to a 
common coordinate system that will be referred to as the 
cluster reference frame. As detailed in Appendix A, the 
transformation for each epoch is derived by minimizing 
the net displacement of a set of "coordinate reference" 
stars, allowing for proper motions, relative to their posi- 
tions in a common reference image, which, in this case, is 
the 2004 July LGSAO image. This procedure attempts 
to ensure that in the cluster reference frame the coordi- 
nate reference stars are at rest (i.e., no net translation, 
rotation, expansion, or skew). A total of '^470 and ~120 
stars serve as coordinate reference stars in the LGSAO 
and speckle epochs, respectively. These stars are selected 
based on the following criteria: (1) high detection cor- 
relations (>0.9), ensuring good positional accuracy, (2) 
located more than O'.'S from Sgr A* to avoid sources with 
measurable non-linear motions (i.e., accelerations in the 
plane of the sky >~8 km/s/yr), (3) low velocities (< 
15 mas/yr, or equivalently ~600 km s~^), which elim- 
inates possible coordinate reference sources that have 



been mismatched across epochs, and (4) lack of spec- 
troscopic identification as a young star from Paumard 
et al. (2006) to eliminate the known net rotation of the 
young stars in the cluster reference frame. Positional un- 
certainties from this transformation process, which are 
characterized by a half sample bootstrap applied to the 
coordinate reference stars, are a factor of ^^1.5 (speckle) 
to 6 (LGSAO) smaller than the centroiding uncertainties 
and grow by less than a factor of 2 between the center of 
the field of view (minimum) and a radius of 3". 

An additional source of positional error originates from 
residual optical distortion in NIRC2. While the residual 
distortion in NIRC2 is small, the extremely precise cen- 
troid measurements in the deep LGSAO images make it a 
significant effect. The presence of such a systematic error 
is established by examining the distribution of positional 
residuals, normalized by measurement (centroiding plus 
alignment) uncertainties, to the linear proper motion fits 
for the coordinate reference stars. The speckle data sets 
do not show large, measurable biases; the speckle mea- 
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Fig. 3. — Comparison of the centroid uncertainties as a function of brightness. Because the very brightest stars (K~9) are saturated 
in their cores in the LGSAO images, there is a slight rise in their centroid uncertainties compared to somewhat fainter sources. Overall, 
however, for bright sources (K<13), the long exposure LGSAO images achieve a centroiding uncertainty of just 0.17 mas, a factor of ~6 
better than the earlier work done with speckle imaging. 



surements, on average, are only la off from the linear 
proper motion fit. In contrast, the much more precise 
deep-LGSAO astrometric measurements are, on average, 
5(7 off from these fits. As described in Appendix B, we 
account for this effect at two stages of our analysis. First, 
0.88 mas is added in quadrature to the positional uncer- 
tainties of the coordinate reference stars to account for 
systematic errors in the coordinate transformations. Sec- 
ond, a local correction, in the coordinate reference frame, 
is derived and applied to the positions of the short period 
stars that were made with LGSAO setups that differ from 
that of the reference image. This procedure ensures that 
residuals from both linear proper motion fits to the co- 
ordinate reference stars (see Appendix A & B) and from 
orbit fits to SO-2 (see §4) are consistent with a normal 
distribution. 

Source confusion can introduce positional biases that 
can be comparable to and, at certain times, larger than 
the statistical errors caused by background or detector 



noise. This occurs when two stars are sufficiently close 
to each other that only one source, rather than two, is 
identified in our analysis with a brightness that includes 
flux from both sources and a position that corresponds 
roughly to the photocenter of the two stars. We divide 
the problem of handling source confusion in our data 
set into the following two cases: (f) the impact of unre- 
solved, underlying stars that are known sources, because 
they were sufficiently well separated at other times, and 
bright enough, to be independently detected, and (2) the 
impact of unresolved, underlying stars that are not iden- 
tified by this study at another time. Because the sources 
are moving so rapidly, instances of the former case are 
easily identified and are typically blended for one year. 
An underlying source that is comparably bright to the 
source of interest can have a significant impact on the 
astrometry; to quantify this effect, we examine the ide- 
alized, noise-free case of a perfectly known PSF by us- 
ing our empirical PSFs to generate idealized binary stars 



7 



and running StarFinder on these simulated images, in- 
putting the known PSF. In this case, the astrometric 
bias is zero once the two components are detected. As 
Figure |4] shows, when the sources are blended, the re- 
sulting astrometric biases can be easily as large as 10 
mas, which is much larger than our centroiding uncer- 
tainties. Such a large astrometric bias occurs when the 
underlying source is at least half as bright as the primary 
source and has a projected, although unresolved, separa- 
tion of ^ 40 mas. We conservatively choose to eliminate 
all astrometric measurements that are known to be the 
blend of two sources from the orbital analysis; specifi- 
cally, if the predicted positions of two known sources are 
separated by less than 60 mas and only one of them is 
detected, then that measurement is removed from our 
analysis. For SO-2 (K=14.2 mag), the eliminated data 
points are those made in 1998, due to confusion with SO- 
19 (K=15.6 mag), in 2002, due to overlap with SgrA*-IR 
(Kmedian = 16.4 mag, but Can be as bright as 14 mag; see 
Do et al. 2008), and in May 2007, due to superposition 
with SO-20 (K=15.9). The impact of these overlapping 
sources, in the first two cases, can be seen in the photo- 
metric measurements (see Figure [5]). 

Source confusion from unknown sources is a smaller 
effect than that from known sources, since the unknown 
stars, in general, are fainter than the known sources. 
Given the long time-baseline of the speckle imaging ex- 
periment, knowledge of sources in this region is most 
likely complete down to K= 16.0 mag. While sources 
as faint as K = 19 mag have been detected in this re- 
gion with LGSAO, crowding and the short time base- 
line of these deeper observations limit the census of 
these sources. Therefore, source confusion from unknown 
sources can give rise to astrometric biases for SO-2 as 
large as 3 mas (from a K=16 mag source), but are typ- 
ically significantly smaller since underlying sources will 
generally be fainter than K = 16 mag. To characterize 
the expected astrometric bias from the undetected source 
distribution, a Monte Carlo simulation was performed by 
generating multiple images with all known stars plus a 
random stellar distribution that, in total, follows the K 
luminosity function and radial profile from Schodel et al. 
(2007). By running these simulated images through our 
data analysis prodedure, we estimate that the astromet- 
ric error from unknown sources for SO-2 is, on average, 0.5 
mas and 1.2 mas for the LGSAO and speckle images, re- 
spectively, and that it scales roughly with the photomet- 
ric bias and galacto-centric distance. However, it should 
be noted that the exact value of this bias is model de- 
pendent. While the photometric bias may be detected in 
the speckle data toward closest approach (see Figure [S]), 
the estimated astrometric biases are smaller than other 
sources of positional uncertainty already included for the 
majority of the SO-2 data points. We therefore do not in- 
corpate them into the reported positional uncertainties. 
Confusion with unknown sources gives rise to larger as- 
trometric biases for SO-16, SO-19, and SO-20, since these 
sources are fainter than SO-2. Given the velocity disper- 
sion in this region and the angular resolution of the data 
sets, the expected timescale associated with biases from 
source confusion is ^^1-2 years. 

As a final step, the relative astrometric positions are 
placed in an absolute coordinate reference frame using 
the positions of seven SiO masers (Reid et al. 2003, 



2007). Infrared observations of these masers with the 
Keck II LGSA0/NIRC2 system between 2005 and 2007 
were obtained with the same camera (i.e., plate scale) 
used for the precision astrometry measurements de- 
scribed above, but with a nine position box pattern and 
a 6" dither offset to create a 22"x22" mosaic of these 
masers (see Appendix C for details). A comparison of 
the maser positions measured in this infrared mosaic to 
the predicted radio positions at this epoch from Reid et 
al. (2003) establishes that the mosaic has an average 
pixel scale of 9.963 ± 0.005 mas/pixel and a position 
angle of north with respect to the NIRC2 columns of 
0.°13 ± 0.°02. This same analysis localizes the radio po- 
sition of Sgr A* in the infrared mosaic to within 5 mas 
in the east-west and north-south directions. By aligning 
the infrared stars detected in both the larger infrared 
mosaic and the precision astrometry image taken during 
the same observing run, we have the necessary coordi- 
nate transformations to convert our relative astrometric 
position measurements into an absolute reference frame. 
For the orbit analysis described in §4, the uncertainties 
in this transformation are applied only after model or- 
bits have been fit to the relative astrometry and are a 
negligible source of uncertainty in the final mass and Ro 
estimates. 

3.2. Spectral Analysis & Radial Velocities 

In the analysis of the spectral data, we accomplish 
the initial basic data processing steps using standard 
IRAF procedures, for NIRC2 and NIRSPEC, and a facil- 
ity IDL data extraction pipeline for OSIRIS. Specifically, 
each data set is first (1) flat fielded, (2) dark subtracted, 
(3) bad pixel and cosmic ray corrected, (4) spatially de- 
warped, and (5) wavelength calibrated. Wavelength cal- 
ibration is performed by identifying OH emission lines 
from sky spectra and fitting a low-order polynomial func- 
tion to the location of the lines. For the NIRSPEC spec- 
tra, neon emission lines from arc lamps provide the wave- 
length calibration. The accuracy of the wavelength cali- 
bration is -9 km s-i or less for NIRC2 and OSIRIS as 
measured by the dispersion of the residuals to the fit. 
Next, the one-dimensional stellar spectra are extracted 
using a spatial window that covers ~ O'.'l for the two di- 
mensional spectral data sets from NIRC2 and NIRSPEC. 
For the three dimensional spectral data set from OSIRIS, 
an extraction box 0'.'14 x 0'.'14 was used. To correct for 
atmospheric telluric absorption features, each spectrum 
is divided by the spectrum of an A- type star. Prior to 
this step, the A-type star's strong intrinsic Br7 feature is 
removed. In the case of the NIRC2 and OSIRIS observa- 
tions, this correction is done with observations of a G2V 
star, which is divided by a model solar spectrum. The 
Br7 corrected region in the G star is then substituted into 
the same region of the A star (Hanson et al. 1996). In 
the case of the NIRSPEC observations, the A-type star's 
Br7 feature is corrected with a model spectrum of Vega^ 
rebinned to the resolution of the A-type star's spectrum 
and convolved with a Gaussian to match the spectral res- 
olution of the observations. The resulting stellar spectra 
are corrected for all telluric absorption features; however, 
they are still contaminated by background emission due 

^ Model taken from the 1993 K urucz Stellar Atmospheres Atlas 

l |f tp : //ftp . stsci . edu/cdbs/cdbs 2/grid/k93models/standards/vega_c95 .fits 
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Fig. 4. — Astrometric bias introduced by an unresolved source in the case of a binary star generated and analyzed with a known PSF. 
Two cases are shown: (solid line) PSF from LGSAO image in 2006 May and (dotted line) PSF from speckle image in 1998 July. The 
contour lines show the amount of bias (in mas) introduced by an underlying source of the indicated flux ratio and separation. Once the 
neighboring source is detected, which happens at separations of ~60 mas, the astrometric bias drops to zero in this idealized case. For 1:1 
binaries, pairs with smaller separations can be resolved. This figure shows that biases well above the positional uncertainties (~1 mas) can 
occur due to underlying sources. 

surement of this line has improved by a factor of 5 with 
improved instrumentation. For the deep LGSAO spec- 
troscopic observations, the radial velocity uncertainties 
for SO-2 are typically -20-25 km s~i. 



to the gas around the Galactic center. The local back- 
ground is estimated and removed by subtracting spectra 
extracted from regions that are '-^O'.'l away. Finally, all 
the spectra within each night of observation are combined 
in an average, weighted by the signal to noise ratio. 

Radial velocity estimates are determined for each spec- 
trum on the basis of the location of the Br7 line. While 
a few of our spectra with broader spectral coverage also 
show a weaker He I triplet at 2.116 /im, we do not in- 
corporate measurements from this line, as it is a blend 
of transitions that can bias the resulting radial veloci- 
ties (see Figure [6|). A Gaussian model is fit to each of 
the Br7 line profiles and the wavelength of the best fit 
peak, is compared to the rest wavelength of Xvacuum = 
2.1661 /im to derive an observed radial velocity. To ob- 
tain radial velocities in the local standard of rest (LSR) 
reference frame, each observed radial velocity is corrected 
for the Earth's rotation, its motion around the Sun, and 
the Sun's peculiar motion with respect to the LSR (U 
— 10 km s~^, radially inwards; Dehnen & Binney 1998). 
Since the LSR is defined as the velocity of an object in 
circular orbit at the radius of the sun, the Sun's peculiar 
motion with respect to the average velocity of stars in 
its vicinity should give the Sun's motion toward the cen- 
ter of the Galaxy. The uncertainties in the final radial 
velocities are obtained from the rms of the fits to the 
line profile measurements from at least three indepen- 
dent subsets of the original data set. Figure [7] shows how 
SO-2's Br7 line has shifted over time and how the mea- 



4. ORBITAL ANALYSIS & RESULTS 

4.1. Point Mass Only Analysis 

To derive the black hole's properties, we assume that 
the stars are responding to the gravitational potential of 
a point mass. In this analysis, the 7 properties of the 
central black hole that are fitted are its mass (M), dis- 
tance (Ro), location on the plane of the sky [Xq, Yq) and 
motion (14, Vy, Vz). In addition to these common free 
parameters, there are the following 6 additional free pa- 
rameters for each star: period (P), eccentricity (e), time 
of periapse passage (To), inclination (i), position angle 
of the ascending node (fi), and the longitude of periapse 
(oj). Using a conjugate gradient minimization routine 
that simultaneously fits the astrometric and radial ve- 
locity measurements, we fit this model to measurements 
that are given in Tables [3] &|11 which includes 27 epochs 
of astrometric measurements and 11 epochs of radial ve- 
locity (RV) measurements, as well as 5 additional epochs 
of radial velocity measurements reported in the litera- 
ture (Eisenhauer et al. 2003, 2005). This excludes all 
the astrometric measurements of SO-2 that are confused 
with another known source (see §3.1). While the 2002 
astrometric data are eliminated due to confusion with 
SgrA*, the 2002 RV points are not, since SgrA* is fea- 
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Fig. 5. — Photometric measurements vs. time for SO-2 (top) and SO-16 (bottom). Measurements that were made when these sources 
coincided with another known source are plotted as unfilled points and excluded from the model fitting procedure. S0-f6 is more affected 
by underlying sources, because it is fainter. Even without a priori knowledge of the underlying sources, their effect is clearly visible in 
photometric measurements made in 1998 and 2002, for SO-2, and 1996-1999 & 2000, for SO-16. 



turclcss and therefore does not bias the measurement of 
RV from SO-2's Br7 absorption Una. In total, there are 
38 astrometric data points and 5 RV measurements. All 
values reported for each parameter are the best fit values 
obtained from minimizing the total x^, which is the sum 
of the from each data type (i.e., xLt = xLt + x|.y); 

The uncertainties on the fitted parameters are esti- 
mated using a Monte Carlo simulation, which is a robust 
approach when performing a fit with many correlated 
parameters. We created 10^ artificial datasets (Nsim) 
containing as many points as the observed dataset (as- 
trometry and radial velocities), in which each point is 
randomly drawn from a Gaussian distribution centered 
on the actual measurement and whose 1 a width is given 
by the associated uncertainty, and run the x^ minimiza- 
tion routine for each realization. Nsim was set to 10^ 
in order to achieve '^6% acctiracy in the resulting esti- 
mates of the a 99.73% confidence limits {3a equivalent for 
a gaussian distribution) of the orbital parameters. Be- 



cause the x^ function contains many local minima, each 
realization of the data is fit 1000 times {Nseed) with dif- 
ferent seeds to find the global minimum. The resulting 
distribution of 10^ values of the fitted parameters from 
the Monte Carlo simulation, once normalized, is the a 
joint probability distribution function of the orbital pa- 
rameters {PDF{0), where O is a vector containing all 
the orbital parameters, Oi). For each orbital parame- 
ter, PDF{0) is marginalized against all other orbital 
parameters to generate a PDF{Oi). The confidence lim- 
its for each parameter are obtained by integrating each 
PDFiPi) from its peak^° outwards to a probability of 
68 

Compared to all other stars at the center of the Milky 
Way, SO-2 dominates our knowledge of the central black 

While the best values from minimizing can differ slightly 
(but well within the uncertainties) from the peak of the PDFiOi) 
values, this has negligible impact on the reported uncertainties. 
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Fig. 6. — The weighted average of all SO-2 spectra obtained with the W. M. Keck II telescope. Since only some of the data sets contain 
the shorter wavelengths, the signal to noise ratio is lower at wavelengths shortward of 2.13 ^m. While Br7 and He I lines are clearly 
detected, only the Br7 line, which is stronger and not the blend of multiple lines, is used to measure the radial velocity of SO-2 as a function 
of time. 



hole's properties. Two facts contribute to this effect. 
Most importantly, it has the shortest known orbital pe- 
riod (P = 15 yr; Schodel et al. 2002, 2003; Ghez et 
al. 2003, 2005a). Furthermore, among the known short- 
period stars, it is the brightest star and therefore the 
least affected by stellar confusion (see Figure 1). Sev- 
eral other stars, in principle, also offer constraints on the 
black hole's properties. In particular, SO-16 is the next 
most kinematically important star, as it is the only other 
star that yields an independent solution for the black 
hole's properties. However, independent solutions for the 
black hole's position from fits to SO-2 and SO-16 mea- 
surements differ by more than 5 a (see Figure [8]). While 
SO-16's measurements in 2000 have already been omitted 
due to overlap with the position of SgrA*, three indepen- 
dent lines of reasoning lead us to believe that some of SO- 
16's remaining astrometric measurements must be signifi- 
cantly biased by radiation from unrecognized, underlying 
stars. First, as shown in Figure IH unknown sources can 
introduce astrometric biases as large as 9 mas for SO-16 
(K=15), in contrast with only 3 mas for SO-2 (K ~ 14), 
because it is only 1 mag above the completeness limit for 
detection in the speckle data set (K'^16 mag; see §3.1). 
Second, a comparison of the solution for the position of 
the black hole (Xq, Yq) based on both the astrometric and 
radial velocity measurements to that based on astrome- 
try alone (fixing the distance, which cannot be solved 
for without radial velocities) yields a consistent position 
from modeling the two cases for SO-2, but produces differ- 
ent results for the two cases from modeling SO-16's mea- 



surements, with the inferred Xq and Yq from astrometry 
alone shifting further away from that obtained from mod- 
eling SO-2's orbit prediction and thereby increasing the 
discrepancy to lOcr. Third and last, while the position 
of the dynamical center from SO-2's orbit is statistically 
consistent with SgrA*-Radio/IR, which is the emissive 
source associated with the central black hole (e.g., Melia 
& Heino 2001; Genzel et al. 2003a; Ghez et al. 2004; 
2005b; Hornstein et al. 2007), the solution from SO-16 
is not (see Figure |8]); this difference cannot be explained 
by allowing the black hole to move with time or by in- 
troducing an extended mass distribution. We therefore 
restrict our remaining analysis to SO-2. 

As shown in Figures iQlfcfTOl the astrometric and radial 
velocity measurements for SO-2 are well fit by a simple 
Keplerian model. For a 13 parameter model (right-hand 
side of figures), the best fit to the data produces a total 
of 54.8 for 57 degrees of freedom (do/) and a /dof 
of 0.961. From the Monte Carlo simulation, we derive 
probability distributions for the central black hole's prop- 
erties, which are shown in Figure[TT]and characterized in 
TableO These distributions give a best fit for the central 
black hole's mass of M^^ = 4.1 ± 0.6 x 10^ Mq and dis- 
tance of Ro — 8.0 ± 0.6 kpc (all quoted uncertainties are 
68% confidence values) . The position of the black hole is 
confined to within ± 1 mas (~100 Schwarzschild Radii). 
As can be seen in Figure [TTl the inferred black hole's 
mass is highly correlated with its distance. Estimates 
from orbital modeling are expected to have a power law 
relationship of the form Mass oc Mq Distance"' with a 
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Fig. 7. — Measurements of SO-2's Br7 line . These three measurements show that, over time, SO-2's radial velocity has changed by 
more than 2,600 km s~^. In order to improve the line detection in the low SNR NIRSPEC observation, the emission from the local gas 
was not removed, which leaves a large Br7 emission feature centered at small radial velocities compared to that of the star at this time. 
With improvements in the adaptive optics system and instrumentation (from NIRSPEC/NGSAO [bottom], to NIRC2/NGSAO [middle] 
and finally to OSIRIS/LGSAO [top]), the precision with which the Br7 absoprtion line can be measured in SO-2 has improved by a factor 
of 5. 

between 1 and 3. For the case of astrometric data only, 
a should be 3 and, for the case of radial velocity data 
only, a is expected to be 1. Currently, the relationship is 
M = (4.1 ±0.1 X 10'^Mq){Ro/8.0 kpcf-^, which suggests 
that the astrometric and radial velocity data sets are 
having roughly equal affect in the model fits for mass"'^^. 

A fit that includes the biased astrometric data points 
significantly alters the best fit solution for SO-2. Includ- 
ing both the 1998 and 2002 data points, which corre- 
spond to confusion with SO-19 and SgrA*-IR respectively, 



The uncertainty in the mass scaling relationship is obtained 
for the case in which Ro is fixed to 8.0 kpc and therefore does not 
include the uncertainty in Rq. 



results in a higher mass (S.TxlO^Af©), distance (9.4 kpc), 
and x^ldof (1.7). Including the 2002 but not the 1998 
data points also produces elevated values (5.2 x IQ^Mq 
and 9.1 kpc) and /dof (1.1). This demonstrates that 
it is important to account for the astrometric biases in- 
troduced by unresolved sources. 

Formal uncertainties in mass and distance estimates 
from orbital fits can be reduced by adding a priori in- 
formation. In particular, it is, in principle, possible to 
constrain the dynamical center to be at the position of 
SgrA*-IR. However, as shown in Figure [51 the six mea- 
surements of SgrA*-IR's position in the deep LGSAO 
images (2005-2007), which have the most precise astro- 
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TABLE 3 

Summary of Keck Astrometric & Photometric Measurements 



UT Date Kobe (mag) X (mas)'' Y (mas)'' 



1995.439 


14.21 




0.09 


42.6 




1.0 


164.10 




0.98 


1996.485 


14.05 


± 


0.10 


53.0 


± 


9.5 


155.4 


± 


9.5 


1997.367 


14.09 




0.07 


56.5 




1.7 


137.0 




1.7 


1999.333 


13.98 


± 


0.10 


66.6 


± 


3.1 


91.5 


± 


3.1 


1999.559 


14.12 


± 


0.04 


67.4 


± 


1.4 


88.3 


± 


1.4 


2000.305 


13.98 


± 


0.14 


64.3 


± 


3.0 


65.8 


± 


3.1 


2000.381 


14.12 


± 


0.04 


66.7 


± 


1.1 


63.0 


± 


1.1 


2000.548 


14.02 


± 


0.08 


64.84 


± 


0.78 


57.94 


± 


0.80 


2000.797 


14.09 


± 


0.11 


65.4 


± 


4.8 


46.8 


± 


4.9 


2001.351 


14.11 


± 


0.08 


56.7 


± 


1.6 


26.5 


± 


1.6 


2001.572 


14.09 


± 


0.04 


53.0 


± 


1.4 


14.2 


± 


1.3 


2003.303 


14.19 


± 


0.04 


-34.9 


± 


1.5 


69.5 


± 


1.6 


2003.554 


14.21 


± 


0.03 


-35.45 


± 


0.90 


81.04 


± 


0.90 


2003.682 


14.28 


± 


0.05 


-34.5 


± 


2.3 


87.4 


± 


2.3 


2004.327 


14.22 


± 


0.05 


-32.15 


± 


0.84 


113.95 


± 


0.86 


2004.564 


14.21 


± 


0.05 


-28.7 


± 


1.4 


121.3 


± 


1.5 


2004.567 


14.21 


± 


0.02 


-28.4 


± 


1.4 


122.9 


± 


1.4 


2004.660 


14.19 


± 


0.04 


-26.8 


± 


1.1 


125.5 


± 


1.1 


2005.312 


14.15 


± 


0.04 


-18.58 


± 


0.88 


142.43 


± 


0.92 


2005.495 


14.28 


± 


0.02 


-18.6 


± 


1.0 ± 1.1 


145.3 


± 


1.0 ± 2.5 


2005.566 


14.18 


± 


0.05 


-15.3 


± 


1.7 


148.9 


± 


1.8 


2005.580 


14.19 


± 


0.01 


-16.9 


± 


0.23 ± 1.0 


146.8 


± 


0.23 ± 1.5 


2006.336 


14.23 


± 


0.01 


-7.97 


± 


0.13 ± 0.77 


159.82 


± 


0.13 ± 0.66 


2006.470 


14.17 


± 


0.01 


-6.01 


± 


0.14 ± 0.77 


161.57 


± 


0.14 ± 0.66 


2006.541 


14.14 


± 


0.01 


-4.89 


± 


0.17 ± 0.77 


162.26 


± 


0.17 ± 0.66 


2007.612 


14.17 


± 


0.01 


6.88 


± 


0.20 ± 0.77 


173.47 


± 


0.20 ± 0.66 



^ X and Y are the relative positions in the EW and NS direction, with increasing 
values to the E and N, respectively. These values are in our absolute coordinate 
system (i.e., relative to SgrA*-Radio; sec Appendix C), but the uncertainties do 
not include the uncertainties in the absolute coordinate system. Measurements 
that are confused with other known sources are not included in this table. ^ Un- 
certainties from residual distortions in NIRC2 relative to the 2004 July reference 
image are reported separately (the second uncertainty term in the table) and 
should be added in quadrature to the other uncertainty terms to obtain the final 
positional uncertainties; since the 2006 - 2007 LGSAO images are all obtained 
with the same set up, positions from these images have correlated residual dis- 
tortion uncertainties. 



TABLE 4 
Summary of Keck Radial Velocity 
Measurements 



UT Date Radial Velocity (km s~^) 
Observed LSR 



2000.487 


1192 ± 100 


1199 ± 100 


2003.438 


-1556 ± 22 


-1550 


± 


22 


2004.474 


-1151 ± 57 


-1143 


± 


57 


2005.410 


-945 ± 16 


-926 


± 


16 


2005.504 


-853 ± 31 


-850 




31 


2006.391 


-715 ± 21 


-692 




21 


2006.461 


-728 ± 17 


-718 


± 


17 


2006.495 


-699 ± 36 


-695 


± 


36 


2006.497 


-717 ± 37 


-713 


± 


26 


2007.385 


-507 ± 50 


-483 


± 


50 


2007.548 


-502 ± 50 


-506 


± 


50 



metric measurements, have an average value that differs 
from the position of the black hole inferred from SO-2's 
orbit by 9.3 mas and a variance of 3 mas, which is a factor 
of 4 larger than expected from the measured positional 
uncertainties (0.7 mas). SgrA*-IR is located where the 
underlying sources are expected to have the highest num- 
ber density and velocity dispersion, which should induce 
time variable positional biases. SgrA*-IR's average K 
magnitude in these deep LGSAO images is 16.4, which is 
comparable to the completeness limit for sources in this 



region (see §3.1) and which is, consequently, potentially 
subject to large astrometric biases (see Figure |4]). We 
therefore suspect that the measured positions of SgrA*- 
IR suffer from astrometric biases from underlying sources 
and do not use its positions to constrain the model fits. 

Another prior, which has been imposed in earlier or- 
bital analyses of SO-2 for Ro (Eisenhauer et al. 2003; 
2005), is on the black hole's motion relative to the mea- 
surements' reference frame. Setting the three dimen- 
sional velocity to zero and fitting a 10 parameter model 
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Fig. 8. — Comparison of estimates of the black hole's location. Colored contours represent the estimates of the dynamical center from 
model fits to kinematic measurements of SO-2 (K = 14.0; blue) and SO-16 (K = 15.0; red). Black contours show the SgrA*-Radio position. 
All contours are plotted at the 68%, 95%, and 99.7% confidence levels (equivalent to 1, 2, and Scr for a Gaussian distribution). The solid 
black points are all the measurements of SgrA*-IR, (K~16) in the maps used for the astrometric analysis. The discrepancy in the black 
hole's location from SO-16's positional measurements appear to be a consequence of biases from unrecognized, underlying stars and thus 
only SO-2's measurements are used to infer the properties of the central black hole. Likewise, the astrometric positions of SgrA*-IR, which 
is even fainter than SO-16, also may be biased (see discussion in §4.1) and are therefore not used to constrain the orbital model used to fit 
SO-2. 



ix^/dof = 1.3; see left-hand side of Figures [i & [TOl) 
yeilds uncertainties in the black hole's properties that 
are a factor of 2 smaller (Rq = 8.0 ± 0.3 kpc and M^h 
= 4.4 ± 0.3 X lO^M©). However this assumption is not 
justified (see, e.g., Salim & Gould 1999; Nikiforov 2008). 
Introducing Vx and Vy (defined such that positive num- 
bers are motions in the E and N directions, respectively) 
into the fit allows the dynamical center to move linearly 
in time in the plane of the sky with respect to the clus- 
ter reference frame. Such an apparent motion can arise 
from either a physical or a data analysis effect. In the 
case of a physical effect, the black hole could be moving 
with respect to the stellar cluster under the gravitational 
influence of a massive companion or the black hole and 
the cluster could be participating in a mutually oppos- 
ing sloshing mode. In the case of a data analysis effect, 



the reference frame could be non-stationary with respect 
to the position of the dynamical center, which might 
arise if there was a systematic problem in our alignment 
of images. Introducing these two parameters therefore 
provides a way of examining possible systematic refer- 
ence frame problems. Fits to a 12 parameter model {Vz 
fixed to zero) to the data have a minimum x^/dof of 
0.95, uncertainties in the black hole's properties that are 
larger than the 10-parameter model, but smaller than 
the 13-parameter model (Rq 8.4 ± 0.4 kpc and M^h = 
4.5 ± 0.4 X IO^Mq), and an estimate for the black hole's 
motion relative to the central stellar cluster of Vx — -0.40 
± 0.25 mas/yr (17 ± 11 km/sec) and Vy = 0.39 ± 0.14 
mas/yr (16 ± 6 km/sec). Since these relative velocities 
are comparable to the constraints on the IR reference 
frame's motion with respect to SgrA*-Radio (i.e., an ab- 
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Fig. 9. — The best fit to the astrometric and radial velocity data, assuming a Keplerian orbital model. The filled points were included in 
the formal fit, while the unfilled points are measurements that are excluded due to source confusion. Uncertainties are plotted on all points, 
except the unfilled/excluded points (here the uncertainties are comparable to the size of the points) for clarity. {Left) To compare with 
what has been done in the past to estimate Ro, we show the fit to the data with a 10 parameter model, which includes the black hole's mass 
(Mi,|j), distance (Ro), and location in the plane of the sky {Xo,Yo) as free parameters, and which fixes the black hole's three dimensional 
velocity {Vx,Vy,Vz) to zero. This results in a ■x^jdoj ~ 1.4. (Right) The data are better reproduced by a 13 parameter model, which 
includes the black hole's mass (M;,;,), distance (Rq), location in the plane of the sky {Xo,Yo), and three dimensional velocity (Vx, Vy, Vz) 
as free parameters, and results in a x^/dof ~ 0.97. Adding these extra free parameters, and in particular Vz, increases the uncertainties 
in the black hole's properties by a factor of two. 



solute reference franie in which the black hole's position 
is known; see Appendix C), it is important to leave Vx 
and Vy as free parameters, even for the case in which one 
assumes that the black hole has no intrinsic motion mo- 
tion with respect to the cluster. Because the black hole 
is so often assumed to be at rest, we report the complete 
solution for the 12 parameter fit (Vz fixed to zero) in 
Tabic El 

As Figure [T^ shows, the black hole's motion along the 
line of sight with respect to our assumed local standard of 
rest (Vz) dominates the uncertainties in Rq in our 13 pa- 
rameter model. Priors on Vz therefore have a signficant 
impact on the resulting uncertainties. Unlike the plane 
of the sky, the reference frame along the line of sight is 
unlikely to have an instrumental systematic drift, since 
each of the spectra are calibrated against OH lines (see 
§3.2). However, it is possible that there is a residual 
radial velocity offset between the LSR and the SO-2 dy- 
namical center. The Sun's peculiar motion with respect 
to the LSR along the line of sight might differ from the as- 
sumed 10 km s~^; that is, the practical realization of the 
LSR is not on a circular orbit around the Galactic cen- 
ter as might occur due the bar potential or to the spiral 



perturbations, so that the average velocity of stars in the 
solar vicinity might have a (small) net radial component. 
Alternatively, the dynamical center of SO-2 could differ 
from the dynamical center of the Galaxy as determined 
at the Sun's (i.e., LSR's) distance, as might result from 
the presence of an intermediate mass black hole compan- 
ion. From the model fit, the implied motion of the LSR 
along the line-of-sight with respect to SO-2's dynamical 
center is -20 ± 33 km/sec, which is consistent with no 
net motion. While no significant motion is detected in 
Vx, Vy, or Vz, the 3cr upper limits for the magnitudes 
of all three are comparable to one another in our 13 pa- 
rameter model (48, 30, and 119 km/sec, respectively). 
Since there are no direct contraints on these quantities 
that can improve these limits, we have allowed them to 
be fully free parameters. However, if we assume that the 
black hole is stationary with respect to the Galaxy, we 
also need to consider the case of Vz set to zero^^. 

Allowing for the uncertainty in the LSR in Vz (it 2 km/sec; 
Gould 2004) produces results that are not distinguishable from 
those reported for the Vz = case. 
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Fig. 10. — The residuals to the best fit Keplerian orbital models shown in Figure 9. The filled points were included in the formal fit, 
while the unfilled points are measurements that are excluded due to source confusion. 



TABLE 5 

Orbital Elements for SO-2 and the Implied Black Hole Properties 



Parameter (Symbol) [Unit] 



Yz =0 Case'' 



Unconstrained Case 



Distance (Ro) [kpc]*" 


8.36 


_l_ 0.30 
=■= 0.44 


7.96 


_l_ 0.57 


Period (P) [yrs] 


15.78 


± 0.35 


15.86 


±8:18 


Semi-major axis (a) [mas] 


124.4 


1 2.4 
=■= 3.3 


126.5 


Eccentricity (e) 


0.8866 ± 0.0059 


0.8904 ± o-oosi 


Time of Closest Approach (To) [yr] 


2002.3358 ± ° 


2002.342 


± 0.010 


Inclination (i) [degrees] 


135.3 


d= 1.3 


134.6 


± 1.3 


Position Angle of the Ascending Node (H) [degrees] 


225.9 


± 1.3 


226.44 


1 0.71 

±}:I 

^ 0.87 


Angle to Periapse (oj) [degrees] 


65.18 


± 1.2 


66.0 


X Dynamical Center (Xq - '>^SgrA,-Radio) [mas]'''"^ 


0.95 


_l_ 0.46 
=■= 1.4 


1.49 


Y Dynamical Center (Yo - Y SgrA,-Radio) [mas]'''"^ 


-4.8 


_i_ 2.2 
=•= 1.6 


-5.4 


± 2.0 


X Velocity [mas/yr] 


-0.40 


± 0.25 


-0.47 


4_ 0.12 
^ 0.33 


Y Velocity (Vy ) [mas/yr] 


0.39 


_i_ 0.09 
=■= 0.18 


0.36 


± 0.12 


Z Velocity (Vi) [km/sec] 






-20 


_|_ 29 


Mass (Mi,fe) [10« Mq] 


4.53 


1 13:1 


4.07 





Density (p) [10^^ Mqpc-3] 


5.83 


^ 0.97 


6.3 


±8:Ii 
_i_ 4:649 

=■= 0.071 


Periapse Distance (Rmin) [mpc] 


0.570 


± 0.037 


0.535 



Note. — Parameters below the horizontal line are derived from those above the line and are provided 
for convenience. 

" Allowing for the uncertainty in the LSR in Vz (± 2 km/sec; Gould 2004) produces results that 
are not distinguishable from those reported above for the Vz = case. ^ The reference time for the 
position of the black hole, when the velocity is a free parameter, is 2000.0'^ Uncertainties in the position 
of SgrA*-Radio are not incorporated into the uncertainties of Xq and Yg. 
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Fig. 11. — The correlation of the estimated black hole's mass and distance. The density of solutions from the MC simulations are shown 
as a color image, with the contours marking the 68%, 95%, and 99.7% confidence limits. While mass and distance are well determined 
from the orbit of SO-2, they are not independent quantities. The exact scaling depends on the relative impact of the astrometric and radial 
data on the model fits. Currently, the inferred mass scales with the inferred distance as M oc Ro^'*. 
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Fig. 12. — Correlation of the estimated black hole's distance and line-of-sight velocity (Vi) from our 13 parameter model fit. Vz 
dominates the uncertainties in Rq and consequently M(,(i. Priors on Vz can reduce the uncertainties in Rq by a factor of two. All previous 
studies have set Vz to zero, which implicitly assumes that there are no massive companions to our Galaxy's central supcrmassivc black hole 
and that the local standard of rest is perfectly known. 



4.2. Point Mass Plus Extended Mass Distribution 
Analysis 

Limits on an extended mass distribution within SO- 
2's orbit are derived by assuming that the gravitational 
potential consists of a point mass and an extended mass 
distribution, and allowing for a Newtonian precession of 
the orbits (see, e.g., Rubilar & Eckart 2001). In order 
to do this, we use the orbit fitting procedure described 
in Weinberg et al. (2005), and adopt an extended mass 
distribution that has a power-law density profile p{r) = 
Poif/ro) This introduces two additional parameters 
to the model: the normalization of the profile and its 
slope 7. The total enclosed mass is then given by 

M (< r) = A/bh + Mcxt(< '^o) f ^ j , (1) 

where we quote values for the normalization Mext(< tq) 
at ro — 0.01 pc, corresponding to the characteristic scale 
of the orbit. Figure fT3l shows the constraint on Afcxt(< 



0.01 pc) and 7 from a fit to the astrometric and radial 
velocity measurements for SO-2. The 99.7% confidence 
upper-bound on the extended mass is Mext(< 0.01 pc) ~ 
3 — 4 X IO^A/q and has only a weak dependence on 7. 

Mouawad et al. (2005) report a similar upper-bound 
on the extended mass in fits to the orbit of SO-2. Their 
analysis differs only slightly from that presented here in 
that it forces the focus to be at the inferred radio posi- 
tion of Sgr A*, assumes a Plummer model mass distri- 
bution, and is based on data presented in Eisenhauer et 
al. (2003). Similarly, Zakharov et al. (2007) use an or- 
der of magnitude analysis to show that if the total mass 
of the extended matter enclosed within the SO-2 orbit 
is > 10^ Mq, then it would produce a detectable apoc- 
enter shift A0 > 10 mas (see also § 3.2 in Weinberg et 
al. 2005). Hall & Gondolo (2006) fit the total measured 
mass concentration M{< r) given in Ghez et al. (2005) 
assuming a power-law density profile and obtain an up- 
per bound of w IO^Mq between 0.001 — 1 pc. 
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Fig. 13. Limits on the amount of mass in an extended distribution contained within SO-2's apoapse distance. The three hnes correspond 
to the 68.3%, 9-5.4%, and 99.7% upper-bound confidence Hmits. The 99.7% confidence upper-bound of 3-4 X 10'' Mq is a fairly weal< function 
of the slope of the assumed power-law mass profile. Simple models of the stellar distribution suggest Mext(< O.Olpc) ~ W^Mq, a factor 
of ~100 smaller than the current measurement uncertainty. 

The surface brightness of stars as a function of pro- 
jected radius from Sgr A* is well measured down to a 
radius of ~0"5 (~0.02 pc). With an assumed constant 
mass to light ratio, the inferred stellar mass distribution 



between this inner radius and an outer radius of 10" is 
consistent with 



M,(< r) = (6 X IO^Mq) (^^^ 



1.6 



(2) 



(Genzel et al. 2003b; see also Schodel et al. 2007). 
Extrapolating this profile down to a radius of 0.01 pc 
gives an enclosed mass M,(< 0.01 pc) w 1 — 2 x 10^ Mq. 
Furthermore, theoretical estimates of the density of cold 
dark matter halo particles suggest that ^ lOOOAf0 of 
dark matter might reside in the inner 0.01 pc of the GC 
(Gondolo & Silk 1999; Ullio et al. 2001; Merritt et al. 
2002, Gnedin & Primack 2004). Likewise, the mass con- 
tribution from a cluster of stellar remnants, as predicted 
by Morris (1993) and Miralda-Escude & Gould (2000), is 



expected to be ~ lOOOAf© within 0.01 pc. Unfortunately, 
these estimates are all smaller than the current upper- 
bound by a factor of » 100. Measurements of stellar or- 
bits with a next generation large telescope are, however, 
expected to be sensitive to an extended mass distribution 
of mag nitude M*(< 0.01 pc) w lO^M© (Weinberg et al. 
2005). 

5. DISCUSSION 

Orbit modeling of astrometric and radial velocity mea- 
surements of short period stars provides a direct estimate 

of the Milky Way's central black hole mass and distance. 
Our analysis of SO-2's orbit yeilds a black hole mass of 
Mbh = 4.1 ± 0.6 X lO^^'Afo'and distance of Rq = 8.0 
± 0.6 kpc, if nothing is assumed about the black hole's 
intrinsic motion. If we assume that the black hole has 
no intrinsic motion relative to the central stellar cluster 
(i.e., no massive companion), but still allow for system- 
atics in the reference frames, then we obtain M^h = 4.5 
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± 0.4 X IO^Mq and distance of Rq = 8.4 ± 0.4 kpc. This 
study shows that there are three systematic errors that 
must be accounted for to obtain accuracy in estimates 
of orbital parameters and this leads to larger uncertain- 
ties than have been reported in the past. First, since a 
dominant source of systematic error in the data set ap- 
pears to be source confusion (see §3 & 4), we use only 
data from the brightest short orbital period star, SO- 
2, and only those measurements that are not confused 
with other known sources. Second, the motion of the 
black hole relative to the measurements' reference frame 
should be left as a free parameter, to account for both 
any possible intrinsic motion of the black hole as well as 
systematics in the astrometric or spectroscopic reference 
frames. Third, while SgrA*-IR is detected with a pre- 
cise position in deep LGSAO images, it appears to be 
biased; therefore, the position of the black hole should 
be treated as a free parameter in the fits in spite of the 
temptation to reduce the degrees of freedom with this 
detection. Because these systematics were not incopo- 
rated into earlier simultaneous estimates of Mbh and Rq 
from the orbit of SO-2, the uncertainties in these initial 
studies were significantly underestimated; Eisenhauer et 
al. (2003, 2005), who do not account for the first two 
systematics, obtain M^h — 3.6 ± 0.3 x 10^ Af© and Rq = 
7.6 ± 0.3 kpc. Ghez et al. (2005a) used SO-2, SO-16, 
and SO-19 simultaneously, and allowed T4 and Vy to be 
free parameters, to derive a mass at fixed Ro of 3.7 ± 0.2 
X 10^ (R-o/ 8 kpc)^ M0, which was pulled down by the 
two astrometrically-biased fainter stars, while Ghez et al. 
(2003) obtained a mass estimate of 4.1 ± 0.6 x 10^ (Rq/ 
8 kpc)^ M0 from SO-2 alone. If we ignore the first two 
effects in model fits to our data, as was done by Eisen- 
hauer et al. (2003, 2005; the only other work to estimate 
Rq from orbits), we obtain a poor quality fit {x^/dof 
= 2.0), uncertainties that are a factor of 2 smaller, and 
somewhat higher values than what we report in Table [5] 
(Mbh = 4.7±0.3 X lO^M© and Rq = 8.6 ± 0.2 kpc). The 
removal of biased astrometric points dominates the shift 
in the black hole's mass and distance to lower values in 
our analysis. This is somewhat suprising as this would 
suggest that similar removal of biased points might lower 
the Eisenhauer et al. (2005) results. However the biases 
may differ, as their early astrometric data measurements 
were made at three times lower angular resolution. An 
astrometric reference frame drift could also explain this 
effect, since Vx and Vy were held fixed in their analysis. 
The addition of Vz as a free parameter dominates the 
resulting uncertainties. In summary, in order to get an 
accurate measure of M^bh and Rq from modeling of the 
short period orbits at the Galactic center, it is critical 
to account for the three sources of systematics described 
above. 

The black hole mass measured here from a stellar orbit 
is larger than the ~ 2 — 3 x IO^Mq inferred from using 
projected mass estimators, which rely on measured ve- 
locity dispersions (e.g, Eckart & Genzel 1997; Genzel et 
al 1997; Ghez et al. 1998; Genzel et al. 2000; see also 
Chakrabarty & Saha 2001). This difference most likely 
arises from the assumptions intrinsic to the use of pro- 
jected mass estimators. In particular, the projected mass 
estimators are based on the assumption that the entire 
stellar cluster is measured, which is not the case for the 
early proper motion studies as their fields of view were 



quite small (r ~ 0.1 pc). Such pencil beam measurements 
can lead to significant biases (see discussions in Haller et 
al. 1996; Figer et al. 2003). An additional bias can arise 
if there is a central depression in the stellar distribution, 
such as that suggested by Figer et al. (2003). These 
biases can introduce factors of 2 uncertainties in the val- 
ues of the enclosed mass obtained from projected mass 
estimates and thereby account for the difference between 
the indirect mass estimate from the velocity dispersions 
and the direct mass estimate from the orbital model fit 
to SO-2's kinematic data. 

A higher mass for the central black hole brings our 
Galaxy into better agreement with the Mbh — c relation 
observed for nearby galaxies (e.g., Ferrarese & Mcrritt 
2000; Gebhardt et al. 2000; Tremaine et al. 2002). For 
a bulge velocity dispersion that corresponds to that of 
the Milky Way ('-lOS km s^^; Tremaine et al. 2002), 
the Mbh — c relationship from Tremaine et al. (2002) 
predicts a black hole mass of 9.4 x 10^ Mq, which is 
a factor of 5 larger than the value of the Milky Way's 
black hole mass used by these authors (1.8 x 10^ Mq 
from Chakrabarty & Saha 2001). The black hole mass 
presented here of 4.1 ± 0.6 x IO^Mq brings the Milky 
Way more in line with this relationship. With one of the 
most accurate and lowest central black hole masses, the 
Milky Way is, in principle, an important anchor for the 
Mbh — f relationship. However, the velocity dispersion 
of the Milky Way is much more uncertain than that of 
other nearby galaxies. Therefore our revised mass has 
only modest impact on the coefficients of the Mbh — f 
relation. 

Revision of the central black hole's mass and distance 
can also, in principle, impact our understanding of the 
structure within our galaxy both on small and large 
scales. On the large scale, if we assume that the black 
hole is located at the center of our Galaxy, then its dis- 
tance provides a measure of Rq. Its value from this study 
is consistent with the lAU recommended value of 8.5 kpc 
as well as the value of 8.0 ± 0.5 kpc suggested by Reid 
(1993), based on a "weighted average" of all prior in- 
direct measurements of Rq . Combining the value for Rq 
from this study with the proper motion of Sgr A* along 
the direction of Galactic longitude measured with VLBA 
in the radio quasar reference frame (Reid & Brunthaler 
2004; fisgrA*.iong — -6.379 ± 0.026 mas yr~^) and the 
Sun's deviation from a circular orbit (Cox 2000; 12 km 
s~^) in the direction of Galactic rotation, we obtain an 
estimate of the local rotation speed, 9q, of 229 ± 18 km 
s~^, which is statistically consistent with other measure- 
ments; these include a value of 222 ± 20 km from the 
review of Kerr & Lynden-Bell (1986) and 270 km s""'^ de- 
rived by Mendez et al. (1999) from the absolute proper 
motions of ^30,000 stars in the Southern Proper-Motion 
survey. As two of the fundamental Galactic constants, 
Ro and Oq are critical parameters for axisymmetric mod- 
els of the Milky Way. Under the assumption that the 
stellar and gas kinematics within our Galaxy are well 
measured, the values of Rq and 9o determine the mass 
and shape of the Milky Way (Oiling & Merrifield 2000; 
Oiling & Merrifield 2001). Of particular interest is the 
value of the short-to-long axis ratio of the dark matter 
halo, q, as it offers a valuable opportunity to distinguish 

consensus value with consensus errors 
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between different cosmological models. As Oiling & Mer- 
rifield (200f) demonstrate, the uncertainty in q for the 
Milky Way is dominated by the large uncertainties in 
Ro and ^o- While our uncertainties in Rq are currently 
too large to constrain q, future precision measurements 
of Ro through stellar orbits may be able to do so and 
could thereby possibly distinguigh between various dark 
matter candidates (Oiling & Merrifield 2001). 

Closer to the black hole, knowing its mass and dis- 
tance from the Sun improves our ability to study the 
kinematics of stars within its sphere of influence. Much 
less kinematic information is needed to determine the 
orbital parameters for stars whose motion is dominated 
by the gravitational influence of the central black hole; 
for instance, with only measurements of a star's posi- 
tion, velocity, and accelaration in the plane of the sky 
along with a single line of sight velocity, a complete or- 
bital solution can be derived once the black hole's mass 
and distance arc well contrained. Improved constraints 
on the central black hole's properties and their degenera- 
cies, as presented here, along with improved astrometry, 
has allowed us to derive orbital information for individ- 
ual stars at much larger galacto-ccntric distances. With 
these measurements, in Lu et al. (2006, 2008), we test for 
the existence and properties of the yomig stellar disk(s), 
proposed by Levin & Beloborodov (2003) and Genzel et 
al. (2003b) from a statistical analysis of velocities alone. 
The direct use of individual stellar orbits out beyond 
a radius of 1" reveals only one, relatively thin, disk of 
young stars (Lu et al. 2008). 

On an even smaller scale, the mass and distance of the 
black hole set the magnitude and time-scale for various 
relativistic effects. Given estimated Keplerian orbital el- 
ements for stars at the Galactic center, we expect to be 
able to measure their stellar orbits with sufScient preci- 
sion in upcoming years to detect the Roemer time delay, 
the special relativistic transverse Doppler shift, the gen- 
eral relativistic gravitational red-shift, and the prograde 
motion of periapse (e.g., Weinberg et al. 2005; Zucker & 
Alexander 2007). These effects will most likely be mea- 
sured with SO-2 flrst, as it has the shortest orbital period 
(P=15 yr), is quite eccentric (c=0.89) and, as one of the 
brighter stars {Kso-2 — 14 mag), it can be measured 
with the greatest astrometric and spectroscopic accuracy. 
The radial velocity signatures of the flrst three effects are 
expected to be comparable to each other and will impart 
a ~200 km/s deviation at closest approach (Zucker & 
Alexander 2007) , when the star is predicted to have a line 
of sight velocity of -2500 km/s based on our updated Ke- 
plerian model. This effect is large compared to the radial 
velocity precision (~ 20 km/sec). Likewise, the expected 
apoapse center shift for SO-2, As = R^^^^Jy-i = 0.9 mas 
(see e.g., Weinberg 1972; Weinberg et al. 2005), is an 
order of magnitude larger than our current measurement 
precision {apos 0.1 mas). Improved adaptive optics 
systems on existing telescopes and larger telescopes (see 
Weinberg et al. 2005) will improve the sensitivity to 
the predicted apocenter shift. To put this measurement 
into context with existing tests of general relativity, it is 
useful to note that one of the strongest constraints on 
general relativity to date comes from the Hulse- Taylor 
binary pulsar, PSR 1913-1-16, which has a relativistic pa- 
rameter at periapse, F = Vsch/f periapse, of only 5 x 10~^, 



~ 3 orders of magnitude smaller than that of SO-2 (Tay- 
lor & Wcisberg 1989; Zucker & Alexander 2007). The 
stars at the Galactic center are therefore probing an un- 
explored regime of gravity in terms of the relativistic ob- 
ject's mass scale and compactness. 

Precession from general relativistic effects also influ- 
ences the timcscalc for resonant relaxation processes 
close to the black hole (see, e.g. Ranch & Tremaine 1996; 
Hopman & Alexander 2006). When precession from gen- 
eral relativity dominates over that from the extended 
mass distribution, the resonant relaxation timescale is 
proportional to Al^f^ x (Jlso/J)'^ >^ P, where J and Jlso 
are the orbital angular momenta for the orbit of inter- 
est and at the last stable circular orbit around the black 
hole, respectively, and P is the orbital period. For a 
given semi-major axis and accounting for the linear mass 

dependence of {Jlso/ J)"^, this results in a M^/^ depen- 
dency. Thus the higher black hole mass inferred from 
this study increases the timescale over which the black 
hole's loss cone would be replenished in the regime where 
general relativity dominates. For the regime where the 
extended mass distribution dominates, the resonant re- 

1 /2 

laxation timescale scales only as MJ^ . A higher black 
hole mass also implies a longer period for the innermost 
stable circular orbit. If the central black hole is non- 
spinning, the innermost stable circular orbit has a period 
4 ix'iQfiMe Periodicities on shorter timescales, 
such as the putative quasi-periodic oscillation (QPO) at 
~ 20 min (Genzel. et al. 2003a; Eckart et al. 2006; 
Belanger et al. 2006) have been interpreted as arising 
from the innermost stable circular orbit of a spinning 
black hole. At the present mass, the spin would have to 
be 0.6 of its maximal rate to be consistent with the possi- 
ble periodicity. However, it is important to caution that 
other mechanisms can give rise to such short periodici- 
ties, such as a standing wave pattern recently suggested 
by Tagger & Melia (2006). Furthermore, claims of a 
QPO in SgrA* have been called into question; Do et al. 
(2008) flnd that the near-IR temporal power spectrum 
of SgrA* is statistically consistent with pure red noise, 
such as might be caused by disk instabilities or intermit- 
tent jet fluctuations, and Belanger et al. (in preparation) 
reach a similar conclusion for the X-rays variations. 

6. CONCLUSIONS 

The short orbital period star SO-2 has been inten- 
sively studied astrometrically (1995-2007) and spectro- 
scopically (2000- 2007) with the W. M. Keck 10 meter 
telescopes. Fits of a Keplerian orbit model to these data 
sets, after removing data adversely affected by source 
confusion, result in estimates of the black hole's mass 
and distance of 4.1 ± 0.6 x IO'^Mq and 8.0 ± 0.6 kpc, 
respectively. While the current analysis is dominated 
by 11 years of astrometric measurements that have ~ 
1.2 mas uncertainties, the LGSAO over the last 3 years 
have positional uncertainties that are an order of mag- 
nitude smaller (100-200 /j,as). With higher strclil ratios 
and more sensitivity. LGSAO measurements arc also less 
affected by source confusion; this is especially important 
for the closest approach measurements, which have to 
contend with source confusion from the variable source 
SgrA*-IR. Following SO-2 for another 10 years should re- 
sult in the measurement of the Sun's peculiar motion in 
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the direction of the Galactic center from the orbit of SO-2 
with a precision of a few km and 1% measurement of 
Rq. At this precision, the measurement of i?o is of par- 
ticular interest because it could reduce the uncertainty 
in the cosmic distance ladder. 
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APPENDIX 
CLUSTER REFERENCE FRAME 

All positional measurements from the individual images {X',Y') are transformed with a full first order polynomial 
to a common reference system {X,Y), which we refer to as the cluster reference frame (see Ghez et al. 1998, 2000, 
2005a; Lu et al. 2008). The transformations are derived by minimizing the net displacements, allowing for proper 
motions, of all the coordinate reference stars (see §3.1) relative to their positions in a common reference image (re/), 
which for this study is the 2004 Jul LSGAO image. Specifically, we minimize the following sum over the coordinate 
reference stars (s): 

s 

where AXs,e = Xs,ref + X (ig - Uef) - Xs^e, AYs^e = Ys^ref + Vy^ X (ig - Uef) - ^s^e, and Ws,e = <^AX^^ + (^AY^ ^ , 

and where X and Y are expressed as the following function of the measured positions X' and Y' for for each epoch 
(e) 

Xe = ao, + ai^ X X'^ + a-i, x F/, 

The coefficients for the reference epoch are fixed to ai^^f — ^2^^/ — 1 and ao^^j = 02^^/ — ^o^c/ — ^irc/ = and 
the coefficients for the remaining epochs (ao^, ai^, a2^, bo^, 61^, and 62^) come from the minimization of D. Because 
of degeneracies between coordinate transformations and proper motions of the coordinate reference stars, the net 
displacment is minimized in two steps. First, D is minimized with the proper motions {Vx and Vy) of the coordinate 
reference stars set to zero in order to obtain preliminary transformation coefficients. Using these initial coeffficients, 
we transform all the positional measurements to a common coordinate system and fit a linear motion model to them in 
order to derive a first pass estimate of the proper motions. Second, D is minimized again, using the preliminary proper 
motions and holding them fixed, while the final transformation coefficients are derived. This procedure produces 
proper motions for the coordinate reference stars that have no significant mean motion. We therefore conclude that 
the resulting cluster reference frame is stable and free of significant systematics. 

This procedure is also used to check the stability of the combined effects of the camera systems and the coordinate 
reference stars. By carrying out transformations that allow for only translation, rotation, and a scale change, we 
examine the apparent stability of the camera's pixel scale and angle relative to that recorded in the header. Figure 
[HI shows, on the left-hand side, that the relative pixel scales for the cameras are stable to within 0.053% (rms) over 
the time baseline of this study and that the uncertainty in the angle relative to the header PA is dominated by 
inaccuracies in the header value (most of the jumps correspond to times when the camera is known to have been 
opened for engineering purposes). This also provides a measure of the resampled NIRC pixel scale relative to the 
NIRC2 pixel scale (1.0269 ± 0.0005) and an absolute NIRC pixel scale of 20.46 ± 0.01 mas pix'^ when combined with 
the absolute NIRC2 pixel scale from Appendix C. On the right-hand side, Figure [Ml displays the results of the same 
excerise but using a set of coordinate reference stars that includes the known young stars; the clear systematic trend 
in the relative pixel scales demonstrates the importance of removing this set of stars with known net rotation from the 
coordinate reference star list. 

NIRC2 GEOMETRIC OPTICAL DISTORTIONS 

Relative stellar positions from the deep LGSAO images have accuracies (~0.2 mas) that are 
an order of magnitude smaller than the currently available optical distortion map for NIRC2 
(|http://www2. keck. hawaii.edu/inst/nirc2/preship_testing.pdf). Since LGSA0/NIRC2 data was obtained 
with four dift'erent setups (e.g., centerings and/or position angles on the sky), imperfections in the optical distortion 
corrections can introduce 1-2 mas systematics, if unaccounted for, into the relative positions of SO-2 (and SO-16). We 
therefore introduce two steps into our analysis to correct for this effect. First, we add, in quadrature, an additional 
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Fig. 14. — The plate scale (top) and the position angle (bottom) over time for all data sets aligned using a set of stars that excludes 
(left) and includes (right) the known young stars. The plate scale is relative to the plate scale in the reference epoch of 2004 July LGS. 
The position angle is the absolute position angle offset from the value reported in the NIRC and NIRC2 instrument headers. Once the 
young stars are excluded, the estimated plate scales for NIRC (squares) and NIRC2 (diamonds) are very stable, approximately 0.05% and 
0.03% (rms), respectively, over multiple years. NIRC shows several systematic jumps in the position angle relative to the value reported in 
the image headers, which is most likely a result of instrument or telescope changes. 

0.88 mas to all the LGSAO positional measurements of the coordinate reference stars, such that the proper motions 
and hence coordinate transformations are not biased. The magnitude of this term is derived by finding the value that 
reduced the average offset of these LGSAO points from the linear proper motion fits, which exclude these points, from 
5(T to la. Second, we derive explicit correction terms for the local optical distortions for SO-2 and SO-16 positions in 
each of LGSAO epochs not obtained with the same set up as the reference image (2004 July), using the orbits of 5 
"calibration" stars (SO-3, SO-7, SO-19, SO-26, and SO-27) that are within 0'.'5 of SO-2. These terms are obtained by first 
using only the speckle data, which are distortion calibrated with respect to the reference image (2004 July/LGSAO; 
see Lu et al. 2008), the reference image (taken with setup#l), and the one other LGSAO image taken with the same 
setup as the reference image to solve for the orbits of the 5 calibration stars. For each LGSAO epoch not included 
in these fits, the average offsets of these five stars' aligned measurements from their predicted location is used to 
characterize the residual distortions for that image (relative to the reference image) at the position of SO-2 and SO-16 
and the standard deviation of the offsets provide an estimate of the uncertainties in these values. Setup #3 is the only 
LGSAO observational configuration, other than that used for the reference image, used in multiple epochs. From the 
measurements with setup #3, it can be seen that the rms of their estimated bias terms (0.24 mas) is smaller than the 
uncertainty in each bias term estimated from the rms of the 5 stars (~ 0.67 mas). This suggests that the bias terms 
are relatively static (see also Appendix A) and that their uncertainties are dominated by our uncertainties in the 
stellar orbits (and possible structure in the distortion on scales <0'.'5). We therefore derive an average bias correction 
value and uncertainty for each setup. The final bias terms, which range in value between 1.6 and 2.6 mas, are added 
to the LGSAO positional measurements made with setups #2-4 in the analysis presented in fJH and their uncertainties 
are added in quadrature with the uncertainties associated with centroiding and coordinate transformation; this bias 
term as already been incorporated into the values and uncertainties reported in Table [H Correlations in the bias 
corrections for setup #3 are applied and accounted for in the Monte Carlo simulations described in §4. 

ABSOLUTE ASTROMETRY 

An absolute astrometric reference frame for the Galactic center was established from radio observations of seven SiO 
masers (Reid et al. 2003, 2007). Relative measurements in the infrared were tied to the absolute frame by observing, 
in the infrared, the red-giant stars that are the source of the maser emission (Figure [TSl fc Table IH])- Observations were 
taken in 2005 June, 2006 May, and 2007 August using LGSA0/NIRC2 (see ^ with 10.86 s integrations in the K' 
band, each composed of 60 co-added 0.181 s exposures in order to avoid saturating the bright masers. A nine position 
dither box pattern was used to construct a 22" x 22" mosaic with two exposures at each position for the 2005 mosaic 
and three exposures at each position for the 2006 and 2007 mosaics. The individual frames for each data set were 
cleaned, undistorted, and then registered and mosaicked using the IRAF xregister and drizzle routines. Subset-mosaics 
were also created with only 1 exposure at each position and were used to derive centroiding uncertainties. StarFinder 
was run on the resulting mosaicked images to extract stellar positions and uncertainties from the RMS error of the 
subset- mosaics. Centroiding errors were typically on the order of 1.4 mas. This yields an IR starlist for each epoch 
with positions in NIRC 2 pixel coordinates. 

The radio maser positions were propagated forward using velocities from Reid et al. 2007 to create a radio maser 
starlist at the epoch of each of the above IR mosaics. Uncertainties in these propogated radio positions are, on average. 
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Fig. 15. — Infrared mosaic measuring the positions of the SiO masers. The 7 masers, whose radio positions are well measured by Reid 
et al. 2007 and which used to establish an absolute reference frame, are circled. Dotted lines depict the outline of example LGSAO (green) 
and speckle (blue) images in which the short period stars are measured and placed in the cluster reference frame. Since the masers sparsely 
sample the area of interest, only low order polynomials are used to calibrate the cluster reference frame (i.e., pixel scale, orientation, and 
position of SgrA*-Radio). 



TABLE 6 
Maser Properites 









IRS 9 


IRS 7 


IRS 12N 


IRS 28 


IRS lOEE 


IRS 15NE 


IRS 17 


Average 


K Magnitude 


9.1 


7.7 


9.5 


9.3 


11.3 


10.2 


8.9 




X Position (arcsec) 


5.679 


0.032 


-3.264 


10.484 


7.684 


1.209 


13.139 




Y Position (arcsec) 


-6.332 


5.529 


-6.912 


-5.833 


4.196 


11.268 


5.560 




X Velocity (mas/yr) 


3.06 


-0.58 


-1.06 


2.00 


0.04 


-1.96 


-1.61 




Y Velocity (mas/yr) 


2.11 


-3.52 


-2.70 


-5.29 


-2.09 


-5.68 


-0.75 






IR - Radio 


X Position (mas) 


3.7 ± 5.2 


3.5 ± 7.2 


11.6 ± 7.3 


0.6 ± 5.9 


-3.9 ± 3.4 


-1.9 ± 5.5 


-5.3 ± 5.2 


1.2 ± 5.7 




IR - Radio 


Y Position (mas) 


-4.1 ± 4.8 


-3.4 ± 7.4 


4.4 ± 5.4 


-4.5 ± 5.9 


0.9 ± 5.0 


9.4 ± 7.5 


-6.4 ± 6.7 


-0.5 ± 5.7 




IR - Radio 


X Velocity (mas/yr) 


1.4 ± 0.8 


0.8 ± 0.5 


1.7 ± 1.0 


-0.3 ± 2.6 


-0.1 ± 0.3 


-1.0 ± 0.4 


1.7 ± 1.2 


0.6 ± 1.1 




IR - Radio 


Y Velocity (mas/yr) 


0.4 ± 0.8 


-1.7 ± 0.6 


0.7 ± 1.0 


-0.5 ± 2.6 


-1.5 ± 0.3 


0.0 ± 0.4 


-3.5 ± 1.3 


-0.9 ± 1.5 



Note. — Mascr positional uncertainties and differences are averaged over the 3 maser epoehs: 2005.495, 2006.336, 2007.612. The values in the last column are the 
average and standard deviation of the values for the individual masers. 
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TABLE 7 
NIRC2 Absolute Astrometry 



2005.495 2006.336 2007.612 

Plate Scale (mas/pixel) 9.963 ± 0.005 9.964 ± 0.004 9.961 ± 0.006 
Angle (deg) 0.12 ± 0.02 0.13 ± 0.02 0.14 ± 0.02 



Note. — The final values for the NIRC2 plate scale, the position angle of 
the NIRC2 columns with respect to North Sgr A* are averages over these values 
and their errors. 

~1.4 mas. For each epoch, the IR maser starhst was ahgned to the Radio mosaic starhst, which resulted in a new IR 
mosaic starhst in the absolute astrometric reference frame with Sgr A*-Radio at the origin. This alignment process 
used only four independent parameters (a global pixel scale, a rotation, and an origin in the x and y directions) to 
transform between the NIRC2 coordinate system of the IR mosaics to the absolute coordinate system of the radio 
masers. While using higher order polynomial transformations reduce the residual offsets positions from SgrA* between 
the infrared and radio measurements, we conservatively chose to use this low order transformation to capture within 
the uncertainties the possible impact of systematics, such as uncorrected residual camera distortions and differential 
atmospheric refraction. This is particularly important given the sparse sampling of masers across the rgion of interest 
(see Figure [T5|) . Uncertainties in the transformation to absolute coordinates, which were determined with a half-sample 
bootstrap Monte Carlo simulation of 100 iterations where each iteration uses only half the stars in each starlist, were 
added in quadrature to the infrared centroiding uncertainties to produce a final uncertainites in the infrared absolute 
positions of the masers. After the transformation to absolute coordinates, the absolute value of the offsets between 
the positions of the masers relative to SgrA* measured in the infrared and radio are on average 0.8cr and 0.8 cr, or 
equivalenly, 5.7 mas and 5.7 mas in the x and y direction, respectively (see TablelB]); we take this to be our uncertainty 
in the position of Sgr A*-Radio in the infrared maser mosaic. Likewise, the transformations between the infrared and 
radio reference frame yields a plate scale of 9.963 ± 0.005 mas/pix and a position angle offset for NIRC2 of 0.13°± 
0.02°. Each of the three infrared maser mosaics yields comparable results (see Table [7]). Uncertainties in the absolute 
positions in the infrared reference frame are dominated by residual optical distortions, which are amplified by the large 
dithers necessary to construct the mosaics. 

A comparison of the maser's proper motions as measured in the radio and the infrared provides an estimate of how 
accurately we can transform our relative measurements into a reference frame in which SgrA*-Radio is at rest and the 
orientation is set by background quasars (Reid et al. 2007). The absolute infrared proper motions of the masers, as 
well as all other stars detected in the infrared maser mosaics, were derived by fitting a linear model to the positions as 
a function of time from the three IR maser starlists that were separately aligned to the radio reference frame. Because 
the alignment uncertainties are dominated by residual distortion and therefore correlated across epochs for a given 
maser, this source of uncertainty is not included in the linear proper motion modeling. The differences in the proper 
motions measured in the radio and in the infrared have an average value of 0.6 ± 0.4 mas/yr and -0.9 ± 0.6 mas/yr 
in the x and y directions, respectively, where the uncertainties are the standard deviation of the mean. Therefore, at 
present, it is not possible to use these measurements to eliminate possible drifts in the cluster reference frame as the 
source of any apparent T4 or Vy from the orbital fits of SO-2 (see §4.1). 

The relative astrometry measurements presented in §3 were transformed into this absolute reference frame through 
a set of infrared stars we designated as infrared absolute astrometric standards. Absolute astrometric standards were 
defined to be those stars that are 1) detected in all three IR mosaics (2005, 2006, 2007), 2) outside the central arcsecond 
(r>0'.'5), 3) have velocities less than 15 mas/yr and velocity errors less than 5 mas/yr, 4) have reasonable velocity 
fits {x^/dof < 4), and 5) are brighter than K=15. With absolute kinematics for 158 stars within 5", we solve for 
a 4 parameter transformation model by comparing the relative positions in the reference epoch image, which are in 
instrumental pixel coordinates, and the estimated absolute coordinates for that epoch, which are in arcsec relative to 
the position of SgrA*. Since all other epochs are aligned to this reference epoch, positional measurements for all stars 
in all epochs are easily transformed into absolute coordinates. While uncertainties in the absolute infrared reference 
frame dominate the final absolute positional uncertainties relative to SgrA*-Radio, they are a negligible source of 
uncertainty for the orbital analysis. 
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